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We study the dynamics of two-component Bose-Einstein condensates trapped in different vertical 
positions in the presence of an oscillating magnetic field. It is shown here how tuning appropriately 
the oscillation frequency of the magnetic field leads to the levitation of the system against gravity. 
This phenomenon is a manifestation of a macroscopic non-adiabatic tunneling in a system with 
internal degrees of freedom. 
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The now classical observation of Bose-Einstein conden- 
sation with dilute atomic vapors in a series of experiments 
has been the stimulus of a great number of experimen- 
tal and theoretical works in the field of matter waves and 
macroscopic quantum dynamics. 

Shortly after the production of single-component Bose- 
Einstein condensates (BECs), multicomponent conden- 
sates were observed 0, in mixtures of two hyperfine 
states of 87 Rb, \F = l,m/ = -1) and \F = 2,m f = 1) 
with and without "sag" corresponding to equal or dif- 
ferent trap centers, respectively. These works, together 
with other experiments which followed soon |j,|5|,|£| stim- 
ulated many studies on multi-component BECs which 
open scenarios different from those found in single com- 
ponent BECs, for instance for their ground state and 
excitations. 

In spinor BECs it is easy to induce Rabi-type tran- 
sitions optically. Similar phenomena involving coupling 
of different states occur in tunneling of BECs in optical 
lattices or double- well potentials [a 0, @, El EH EH E3 

mm. 

In this paper we describe a striking levitation phe- 
nomenon which occurs when both the spatial aspects 
of the dynamics are put together with the two inter- 
nal degrees of freedom present in a two-component Bose- 
Einstein condensate in which transitions between compo- 
nents are driven by an oscillating magnetic field. In this 
scenario the dynamics splits into a fast complex spatio- 
temporal oscillation of the condensate wavefunctions to- 
gether with a slow dynamical levitation of the total center 
of mass against gravity with neither applied mechanical 
force nor associated classical trajectories. 

In this paper we will study the dynamics of two- 
component BECs where each component is trapped in 
different vertical positions as schematically indicated in 




initial location 
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FIG. f: Schematic plot of the setup studied in this paper. 
Each of the components of a two-component BEC is subject 
to the action of gravity and different confining potentials Vi 
and Vi with minima z\ and z%, respectively. An oscillating 
magnetic field induces transitions between both states. A 
Franck-Condon type vertical transition is indicated by the 
black arrow. Initially all atoms are in the first component. 



Fig. ^ An optical dipole trap together with a mag- 
netic field gradient can provide such kind of confinement 
To fix ideas we will consider our system to be the 
\F = 1,771/ = —1) and \F — 2,mf = 1) hyperfine states 
of 87 Rb coupled by an oscillating magnetic field of fre- 
quency £1 inducing transitions between them. We shall 
pay a special attention to the situation of the Franck- 
Condon- type vertical transition 0, 0] shown in Fig.l. 
To simplify the treatment yet preserving the spatial as- 
pects of the dynamics we will assume our system to be 
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magnetically tightly confined along one of the transverse 
directions |16( to two effective dimensions. 
The Hamiltonian density of our system is 



3=1 V 

+ u 12 \^i\ 2 \H 2 
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where Vj are the confining potentials acting on each of 
the species given by 

Vj(x, z) = ^muj 2 [x 2 + (z- z 3 ) 2 ] + ej + mgz, (2) 

where ej stands for the internal electronic energy and 
g is the gravity constant. Vj(x,z) has the minimum 

value mzjg - y(^) 2 + e j at ( x > z ) = (0> 2 j = z j ~ Q/^ 2 ) ■ 
U\\ 1 U22 and [7 12 = U21 are the effective nonlinear in- 
teraction coefficients defined by f/y = Anti 2 ''a^N '/ 'mi in 
2-dimcnsional (2D) space where N is the total number 
of atoms and ay are the scattering lengths for binary 
collisions. Finally, I = y^h/mu: is a characteristic length. 

The Gross-Pitaevskii equation (GPE) for the mean 
field dynamics of our system is derived through Lagrange 
equations with use of the Hamiltonian H = J J Tidxdz. 
Additionally, we change to new variables scaled as u)t — > 
t, x/l — > x, zj I — > z, £ip — > ip and obtain 
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where Vj = Vj(x, z) = \[x 2 + (z - zj) 2 } + g'z + with 
\]'.. = iirNaij/l,B' = B/hj,fl' = Q/Hw,g' = g/£uj 2 and 
e'j = ej/fkv being dimensionless constants. The normal- 
ization condition for tp is J/(|i/'i| 2 + IV^I 2 ) dxdz = 1. 

Initially we prepare a circularly-symmetric Gaussian 
wavepacket with its center of mass at ri = (0, z^), 
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x 2 + (z - z{f 
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which minimizes the energy E = J d 2 r( | |V-0i| 2 + 

^1 l^i 1 2 + ^2 i l' ! /'i| 4 ) over a family of gaussian functions 
and approximates the ground state when all atoms are in 
state |1). We then apply an external microwave field that 
induces a Franck-Condon type transition between the 
species with frequency tt' = AV{ 2 = F 2 '(0, z^-F/fO, z{). 
Without loss of generality, we will consider the case 
Z2 = — z* . After the oscillating magnetic field is switched 
on, we study the dynamics according to Eq. 




FIG. 2: [Color online] Time dependence of the population 
mixing for the case of resonance of Franck-Condon type. 
Dashed and solid lines stand for Pi = //l^il 2 an d Pi = 
// |V>2| 2 , respectively for B' = 50. (a) Sz = 1.0 (b) Sz = 5.0. 
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FIG. 3: (a) Time evolution of the center of mass z = 
/ i/' t zipd 2 r = J z (|^i | 2 + |i/> 2 | 2 ) d 2 r for B' = 50 and Sz = 1.0 
(broken line) and Sz = 5.0 (solid line). The dotted line 
stands for the case far from the standard resonance (fi' = 
l,8e — 500) with Sz — 5. (b) and (c) Pseudocolor plots 
of *(x,z,t)| 2 = \ipi(x,z,t)\ 2 + \4> 2 (x, z, t) I 2 for t = and 
t = 2.08 on the region (x,z) G [-6.5,6.5] x [-9.6,9.6] show- 
ing the rising (white broken lines) of the total wavefunction. 



Since our results are based on long-time numerical sim- 
ulations of Eqs. we have integrated them by differ- 
ent numerical methods: an alternating direction implicit 
method, a Crank-Nicholson's method, and a split-step 
method with pseudospectral integration of the spatial 
derivatives 17] , all of which lead to identical results. The 
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4: Pseudocolor plots of \tp2(x, z,t)\ 2 on the region (x,z) £ [—6.5,6.5] x [—7.1,6] for different values of t and Sz — 5.0, 
50. The time spacing between all pannels, except for pannels 2 and 3 and pannels 5 and 6, is half the Rabi period trb- 




FIG. 5: Surface plots of \4>i(x,z,t = 2.08) | 2 (red) and 
\ip 2 (x,z,t = 2. 08)\ 2 ( green) for the situation described in Fig. 
|1| The arrow indicates the location of maximum density of 
| 2 corresponding to a low-density region of \ip2 \ 2 ■ The spa- 
tial region shown in the plot is (x, z) £ [—6.5, 6.5] x [—6.5, 6.5] 



values of the scaled parameters used throughout the text 
are: U' lx = 100, U' 22 = 97, U{ 2 = 94; g' = 0.1; z\ = -2.5 
and z 2 = ~ z \ with Sz — z 2 — z\ — 1, 3 and 5. e[ — —250, 
e' 2 = 250 and Se' = e' 2 - e[ = 500 []]|. 

First we study the population mixing between the com- 
ponents. In the standard resonance case (Q 1 — 5e' = e' 2 — 
e'i) without coupling with the orbital degrees of freedom, 
the Rabi formula implies that each Pj(t)(= f \ipj\ 2 d 2 r) 
shows a fast oscillation whose frequency is proportional 
to the magnetic field amplitude B' In our Franck- 

Condon type transition (Q' = AV{ 2 ), the spatial orbital 
functions of BEC are coupled with the electronic de- 
grees of freedom, and the Rabi oscillation (with period 
trb = 2tt/B') becomes modulated sooner or later. In 
Fig. [2] we show the time dependence of both populations 
for different values of the distance between trap centers 
5z. We find that unless Sz is sufficiently small, the regu- 
lar oscillation responsible for population mixing becomes 
suppressed as time elapses, which is attributed to the 
energy conversion from electronic to mechanical energy. 

We have also studied the motion of the total center 
of mass (CM) of the system along the vertical direction 
defined by z = J xjj^ zxpd 2 r — J z (\ipi\ 2 + IV^I 2 ) cPr. 

In Fig. |3fa), except for the lowest curve, we show 
typical examples of the time evolution of the center of 
mass. In the Franck-Condon resonance situation de- 



scribed above we find a rising of the total center of mass 
against gravity. This is an interesting phenomenon which 
happens while the underlying population mixing shows 
a fast Rabi oscillation. In our numerical simulations 
we observe that the CM of the condensate exhibits a 
smooth and upward acceleration, reaching a maximum 
height ~ z 2 for t = 2.2 corresponding to (2.2/trb ~ 18 
Rabi periods) for B' — 50. This physical phenomenon 
is a levitation of the condensate through a macroscopic 
quantum transition which is not accompanied by neither 
an applied mechanical force nor a classical trajectory. 
Figs. GUb) and (c) are the total probablility density 
|*(a;, z, t)\ 2 = |V>i(:c, z, t)\ 2 + \ip 2 (x, z, t)\ 2 before and after 
the levitation, respectively. The levitation effect becomes 
more pronounced as Sz is increased. It is very interesting 
to point our that no levitation occurs in the case far from 
the standard resonance (£1' <ti Se'), which is confirmed in 
the lowest curve in Fig. 01a). 

Fig. 0]shows the spatiotemporal dynamics of the wave- 
function \ip 2 (x, z, t)\ 2 for times located in the three re- 
gions which are clearly discriminated in the dynamics for 
Sz = 5. First, in the early stage region (t = ~ 0.5) 
the wavefunction oscillates near the minimum (zl) of the 
lower potential well. As time goes on (t = 1 ~ 1.5), 
the central position of the solutions moves towards the 
minimum of the upper well (z 2 ). Finally, when the CM 
position reaches its peak value ~ z 2 for t = 2 ~ 2.2, 
the wavefunction has a double-humped structure with 
its valley located near the minimum of the upper well. 
Fig. [S]shows that the two components l^i^l 2 form a do- 
main structure in this time region, a typical feature of the 
system due to the nonlinear interactions. The nonlinear- 
ity plays an important role in the formation of domain 
structures during the transitions, a feature which is not 
present in the noninteracting case. 

The origin of this transport phenomenon is that the 
condensate is effectively driven by the oscillating mag- 
netic field to the higher energy states of the upper level, 
thus leading to a Franck-Condon type vertical transition. 
After the transition the new position is unstable, and the 
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atoms relax towards the new minimum (at the position 
with larger z), leading to the observed transport. 

We have studied the dependence of the levitation phe- 
nomenon on the problem parameters B' and Sz = z\ — z\. 
In Fig. IHIa) we show the dependence on B' of the upward 
acceleration rate, a, evaluated by using a quadratic func- 
tion of t (:z — z* + \ol1?) which is fitted to the CM data 
in the early stage region. Interestingly, a has a maximum 
at some optimal value of B' . This fact indicates that the 
Rabi oscillation is essential for the energy conversion be- 
tween electronic and orbital degrees of freedom; however, 
when the Rabi oscillation is too fast to guarantee the re- 
laxation towards z%. there is a reduction of the effective 
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FIG. 6: Upward acceleration rate (a) for the motion of the 
center of mass during the levitation process (a) B' dependence 
in the cases Sz = l(dotted line), Sz — 3 (dashed line), 5z = 5 
(solid line) ; (b) Dependence of a on Sz for B' — 50 the Franck- 
Condon (solid line) and standard (dotted line) resonances. 

In Fig. Etb) we see that under a fixed B', a increases 
monotonically in Sz, but indicates a slight derivation 
from the linear law in the case of non-Franck-Condon 
type transition (e. g. in the case of standard resonance, 
Of = 5e'). 

In conclusion, wc have studied the dynamics of two- 
component Bose-Einstein condensates subject to an a.c. 
driving magnetic field. In the case of Franck-Condon's 
vertical transition in which the driving frequency f2 cor- 
responds to the energy difference between the trapping 
potentials V\ and V2 at the minimum position z* of V\, 
the condensate initially located around z\ is effectively 
driven to unstable states of V2 , relaxing towards the po- 
tential minimum z^ (= z* +Sz) of V2 ■ In this situation the 
population mixing shows a fast Rabi oscillation, while the 
mean center-of-mass of the condensate exhibits a smooth 
and upward acceleration with its rate (a) showing a max- 
imum at some optimal magnetic-field amplitude and a 
monotonic increase in Sz. This phenomenon is a mani- 
festation of macroscopic non-adiabatic quantum tunnel- 
ing with internal degrees of freedom, which differs es- 
sentially from other quantum tunneling processes with 
Bose-Einstein condensates discussed up to now. 
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